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We study optical properties of arrays of ultrathin nanowires by means of the Brillouin scattering 
of light on magnons. We employ the Stokes/anti-Stokes scattering asymmetry to probe the circular 
polarization of a local electric field induced inside nanowires by linearly polarized light waves. We 
observe the anomalous polarization conversion of the opposite sign than that in a bulk medium or 
thick nanowires with a great enhancement of the degree of circular polarization attributed to an 
unconventional refraction in the nanowire medium. 


The study of magneto-optical response of tailored 
nanostructures is in the focus of active research of 
nanostructured materials m- Nonmagnetic metallic 
nanowires are well known in optics, and they are em¬ 
ployed as building blocks of the so-called wire metama¬ 
terials [2]. Such structures demonstrate many unusual 
properties, including negative refraction [S], enhanced 
sensing [6], super-lensing [7], strong nonlocal effects [8], 
nonlinearities [S], and they can boost light-matter inter¬ 
action in the regime of the hyperbolic dispersion m- 

Here, we study magneto-optical (MO) properties of ul¬ 
trathin ferromagnetic Co nanowire arrays (see Fig. Sby 
means of polarization-resolved Brillouin light scattering 
on magnons [IIHI3]. In a sharp contrast to the previous 
studies of ferromagnetic nanowires with the diameter D 
of 20 nm here we study arrays of thin nanowires 
with diameter D ^ 4.8 nm and lattice spacing ~ 17 nm 
(Fig. §. We reveal a striking anomaly of Stokes/anti- 
Stokes pattern in the spectra of Brillouin Light Scatter¬ 
ing (BLS) from magnons in such structures. First, its 
asymmetry is inverted, which is impossible in principle 
in continuous non-structured metal layers; secondly, it 
is unconventionally large. Although the latter effect can 
be predicted formally for continuous films, such a huge 
asymmetry can only be calculated for materials with un¬ 
realistic optical parameters. This is explained by a strong 
modification of optical and magneto-optical properties 
for thin nanowires. While the coherent propagation of 
photons in the sample is weakly affected by thin wires 
due to their small volume concentration, the photon in¬ 
teraction with magnons is confined to the volume inside 
the wires where the wave polarization is strongly modi¬ 
fied due to a huge mismatch in the dielectric constants of 
metallic nanowires and dielectric media. As a result, we 
observe that a linearly polarized light obliquely incident 
upon a metacrystal composed of ultrathin Co nanowires 
acquires ellipticity inside the wires that is enhanced by 


an order of magnitude being of the opposite sign as that 
in the case of continuous Co films or thicker Co wires. 
We visualize this effect by measuring the asymmetry of 
Stokes and anti-Stokes peaks in the Brillouin scattering 
spectra of light by the spin-wave modes of the wires, since 
the scattering can probe local electric fields m- We ex¬ 
pect that our results will be instrumental for the emerg¬ 
ing field of nonlinear spectroscopy of metamaterials [16] 
as well as for a design of novel structures with strong 
chiral responses m and polarization-sensitive light root¬ 
ing [iHHini- 

An array of ultrathin Co nanowires has been grown by 
sequential pulsed laser deposition of Co and Ce02 in re¬ 
ductive conditions (Hgrowth = 10“® mbar), leading to self- 
assembly of Co nanowires embedded in a Ce02 matrix 
on a SrTiO3(001) (SurfaceNet GmbH) substrate using 
a quadrupled Nd:YAG laser (wavelength 266 nm) oper¬ 
ating at 10 Hz and a fluence in the 1-3 Jxcrn”^ range. 



FIG. 1: Schematics of the Co nanowires and the Brillouin 
scattering geometry. Incident (p-polarized) and scattered (s- 
polarized) photons and the direction of the magnetic field 
H are indicated. Left inset: TEM image, acquired in energy- 
filtered mode at the Co L edge, taken in plane-view geometry. 
Right inset: HRTEM image of a single nanowire taken in 
cross-section. 
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More details are given in Refs. mi 121]. Metallic nanowire 
formation in the sample was evidenced using high resolu¬ 
tion and energy-filtered transmission electron microscopy 
data (acquired using a JEOL JEM 2100F equipped with 
a field-emission gun operated at 200 kV and a Gatan GIF 
spectrometer), see the right inset of Fig. The wires are 
oriented perpendicular to the surface of the substrate and 
extend throughout the matrix thickness, t. Importantly, 
the technology employed is such that the wire length h 
coincides with the film thickness t, i.e. t = h. Their 
length turned out to be equal to 470 nm. The diameter, 
D = 2R = 4.8 ± 0.7 nm of the wires was determined by 
collecting images in a planar geometry. 

Now we proceed with the analysis of the spectra of 
Brillouin Light Scattering (BLS) from thermal magnons 
localized on ferromagnetic nanowires. The experimental 
arrangement is sketched in Fig. and it corresponds to 
the Damon-Eshbach geometry [23|. A magnetic field H is 
applied in the plane of the sample. The plane of incidence 
is perpendicular to the applied field. The incident laser 
beam is p-polarized with wavelength Aopt = 532 nm. The 
backscattered light is probed in s-polarization through 
a tandem Fabry Perot interferometer (JR Sandercock 
product). The BLS process is a particular case of non¬ 
linear wave-mixing, and it generates at the output two 
frequency shifted optical waves, namely a down-shifted, 
known both in Raman and Brillouin spectroscopy as the 
Stokes (S) line, and up-shifted called the anti-Stokes (AS) 
line. Typically, light scattering spectra are asymmetric, 
i.e. the amplitudes of S and AS spectral lines are not 
equal. However, the physical mechanisms producing this 
peculiar asymmetry are completely different. In the Ra¬ 
man case it is due to a greater difference between the 
frequencies of S and AS lines which results in an appre¬ 
ciable asymmetry in the density of states corresponding 
to the frequencies w-l-H and w —H, where w and are in¬ 
coming photon and magnon frequencies, respectively. In 
the BLS process the frequency shifts are smaller by sev¬ 
eral orders of magnitude, and entirely different physical 
effects are involved, namely a very particular symmetry 
of MO interactions. Mathematically, the symmetry of 
MO coupling is described by the totally antisymmetric 
Levi-Givita tensor jlllisi. As a result, the efficiency of 
the MO interactions is expressed via the mixed product 
Eg ■ {m xEi) of the polarizations of the interacting waves: 
the incident (Ei) and scattered (Eg) optical waves and 
the scattering spin wave {m). Importantly, in the gen¬ 
eral case of complex vector space this mixed product is 
not invariant with respect to complex conjugation of the 
waves. In physical terms it is the elliptical polarization 
that is linked to the complex vectors and the complex 
conjugation corresponds to the inversion of the direction 
of rotation of such polarization. This crucially important 
point will be revisited below. 

Experimental BLS spectra for the angle of incidence 
9 = 20° are presented in the panel (a) of Fig. The spec- 




Incidence angle 

FIG. 2: (a) Experimental Brillouin scattering spectrum of the 
nanowire array for three incidence angles 9 = 20° (a), with an 
applied in-plane magnetic field of 7 kOe. (b) Spectrum for the 
1.2 nm thick continuous Co film with the field of 6 kOe. (c) 
Measured (squares) and calculated (line) dependence of the 
ratio of Stokes and anti-Stokes spectral peak heights in the 
incidence angle 9. The inset shows the circular polarization 
degree of the electric field inside the wires at j/ = x = 0, 
averaged over the wire length. The calculation details are 
given in text. 


tra are not symmetric: the down-shifted Stokes line and 
the up-shifted anti-Stokes line have different magnitudes, 
which is not untypical of magneto-optical BLS spec¬ 
tra El- What is really not conventional, however, is an 
extremely high degree of the Stokes/anti-Stokes (S/AS) 
asymmetry and even more so its inversion with respect to 
its “classical” pattern in which the intensity of the down¬ 
shifted Stokes peak \Es\'^ is greater than that of the up- 
shifted anti-Stokes one \^s\'^ > {EasI"^- The 

latter is illustrated in Fig. [^b). This dramatic reversal 
of the scattering spectra asymmetry constitutes the main 
result of our work and is analyzed in detail below. 

Now we proceed to the theoretical analysis of the 
Stokes/anti-Stokes asymmetry in the scattering spectra. 
In this respect, identification of the spin wave (SW) 
modes contributing to the BLS spectra is very impor¬ 
tant. The applied static magnetic field of 7 kOe fully 
saturates the sample so that magnetization in each wire 
is perfectly homogeneous. The latter allows a reliable 
theoretical description of the spin-wave behavior, includ¬ 
ing explicit expressions for magnetization. 

The only candidate for the role of the “effectively scat¬ 
tering mode” are the vertical Kittel SW modes, i.e. spin 
waves propagating along the wire axis with the wave 
number K with uniform cross-section distribution (at 
least, in the approximation KR <C I). The limiting case 
of AT = 0 corresponds to a purely magnetostatic perfectly 
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FIG. 3: (a) Theoretically calculated spatial maps of the cir¬ 
cular polarization of the electric field in the unit cell of Co 
nanowire array. The arrows show the propagation directions 
of incident and scattered waves, (b) The cross-sections of the 
spatial map along the wire centers x = y = Q, calculated nu¬ 
merically (black/solid lines) and analytically (red/thin lines), 
(c) Intensity of the electric field averaged over x coordi¬ 
nate within the unit cell a,t y = 0,. The calculation has 
been performed for Aopt = 532 nm, 6 = 20°, TM polariza¬ 
tion, d = 4.8 nm, a — 20 nm, eco = —10.4 -|- 17.li [25] . 
£Ce 02 = 4.8 [26], esrTiOa = 6.0 |27|. Analytical curves corre¬ 
spond to the isotropic Maxwell-Garnett approximation with 
smg = 5 -t- 0.4i. 


uniform oscillations |5| , and it is characterized by the fre¬ 
quency = H{H — 27 rMeff). Here, takes into 

account the dipolar inter-wire interactions between static 
magnetizations in individual nanowires, which makes it 
slightly smaller than the conventional bulk value for the 
cobalt, and 7 is the gyromagnetic ratio. The state of 
ellipticity of magnon polarization is crucial in estimat¬ 
ing the S/AS asymmetry. Its actual value is a result of 
a trade-off between two trends, namely, a purely circu¬ 
lar shape imposed by the ferromagnetic resonance and 
the flattening effect of dipolar interactions. Supplemen¬ 
tal Material includes more details on the dispersion rela¬ 
tion f 2 (Al), spatial distribution and ellipticity of the spin 
waves. 

For theoretical analysis of the optical properties, we 
use the general semiclassical theory of light scatter¬ 
ing |5Sj as sketched below. First, the electric field 
exp(—iwt) of the plane wave incident upon the 
sample at the frequency uj is calculated. The field spatial 
distribution is strongly inhomogeneous and it is modified 
by the interaction with the wires. Second, the electro¬ 
magnetic polarization induced in the wires by the inter¬ 
action with the spin waves is determined [His]: 

Ps,u^-n{r) = , (S) , (1) 

Ps,uj+n{r) = ^mn{r) x , (AS) . ( 2 ) 

Here, ma{r) is the magnetization profile of the given 
spin mode with the frequency H, and ^ is the interaction 
constant. Equations 0 and 0 correspond to Stokes 
and anti-Stokes scattering, respectively. Depending on 


whether the magnon is destructed or created, either m 
or m* enter the expression for polarization. Clearly, the 
absolute values of the polarizations induced at Stokes 
and anti-Stokes processes can differ provided that the 
electromagnetic and spin waves have nonzero ellipticity. 
Finally, the detected held is determined as Es^u}±a{'i') = 
f d^r'G^±n(r,r')Ps^^±a(r') , where G^±n(r,r') is the 
tensor electromagnetic Green function at the correspond¬ 
ing frequency. 

In the experimental geometry, the incident wave is p- 
polarized, i.e. the electric held is in the xz plane (see 
Fig. [^. The scattered wave is detected in s polariza¬ 
tion, i.e. the electric held is parallel to the y axis. 
We have verihed numerically, that \Gyy\ ^ \Gyx\ and 
\Ei^y\ <C \Ei^x\, \Ei,z\ inside the wire regions, i.e. the (lin¬ 
ear) mixing between s and p polarizations inside the wires 
is negligible. Hence, the scattering is determined by the 
y component of the polarizations Eq. 0 and Eq. 0 . 
parallel to the detector polarization. The asymmetry be¬ 
tween Stokes and anti-Stokes scattering can be quantihed 
by the difference of the intensities iHj.w-n.yp— 
that is equal to 


A/s/as = \e\rnn?\Ein?P^i:lyP^Sly , (3) 

where we introduce the coordinate-dependent circular 
polarization degree of the wave e 


p. (r)-—i\e*xe] ^ 2 Ime;(r)e.(r) 


The quantity 0 changes from 1 , for right-circularly po¬ 
larized helds, to — 1 , for left-circularly polarized and it 
vanishes for linear polarization. Equation (|^ demon¬ 
strates, that the Stokes - anti-Stokes asymmetry requires 
the non-zero ellipticity of both interacting waves. For the 
spin wave m polarization is dominated by the ferromag¬ 
netic resonance and is fairly close to circular. Our estima¬ 
tions show that it is of the order of 0.8. Another necessary 
ingredient for the S/AS asymmetry is the non-zero local 

(W''\ 

circular polarization degree of the incident wave Pcircy 
Moreover, to explain the experimentally observed unex¬ 
pectedly strong S/AS asymmetry it is required that it 
is close to the SW ellipticity, i.e. 0.8. To justify the 
observed inversion of the S/AS asymmetry the optical 
ellipticity should be reversed with respect to the case of 
a continuous Co film. Even though the incident electro¬ 
magnetic wave is linearly (p) polarized, the finite ellip¬ 
ticity is induced due to its refraction at the interfaces. 
This effect can be most simply illustrated for the case 
of thick continuous Co film. The p-polarization vector 
inside the sample is equal to Sp oc k^Sx — kx^-z ) where 
kx = wsin 0 /c is the in-plane wave vector determined 
by the incidence angle 9 and k^ = \/{^/cY£co — k"^. 
Since the permittivity of Co at the considered wave¬ 
length A = 532 nm is complex, the z component of 
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the wave vector is complex as well, and the local elec¬ 
tric field is elliptically polarized with Pcirc ~ —0.13 at 
the incidence angle 9 = 20°. A crude Maxwell-Garnett 
model [53] for the Co/Ce 02 nanowire array describes 
it as a slightly lossy dielectric with the averaged per¬ 
mittivity SmG = 2{£xx + £yy)/5 + ^ 5 + 0.4i. 

This Maxwell-Garnett approach yields Pdifjy « —0.012, 
i.e. even smaller and also negative circular polarization. 
Both these predictions are in stark contradiction to ex¬ 
periment. In order to resolve this controversy we have 
resorted to full-wave numerical simulation of the electric 
field profile inside the wires and its circular polarization 
degree using the GST Microwave Studio software pack¬ 
age. The results are presented in Fig. SI panel (a) shows 
the spatial map of the circular polarization degree within 
the array unit cell. Figure [sTj^b) and Fig. [Sl|c) show the 
^-dependence of the polarization at the wire center and of 
averaged electric field (thick/black lines). Thin/red lines 
correspond to the analytical isotropic Maxwell-Garnett 
model. One can see, that the Maxwell-Garnett approx¬ 
imation well describes the distribution of the field and 
the polarization in the air (the region z < 0), governed 
by the interference between incident and specularly re¬ 
flected waves. The field decay in the sample {z > 0) 
along the wire axis is captured by the Maxwell-Garnett 
approximation as well, see Fig. Sl h). However, the cir¬ 
cular polarization inside the wires is strongly different 
from the effective medium model. Contrary to the naive 
effective medium approximation, the circular polariza¬ 
tion inside the wires has a positive sign and is quite 
large. Namely, Pcirc oscillates around the value -1-0.59 
which greatly exceeds the values both for continuous 
Co film (Pcirc ~ —0.13) and for the effective medium 
(-Pcirc = —0.012). These numerical findings fully explain 
our experimental data. 


The polarization reversal can be qualitatively under¬ 
stood in the geometrical optics approximation by exam¬ 
ining the refraction of a plane wave at interfaces of the 
structure, as shown in Fig. |^a,b). In the case of con¬ 
ventional positive refraction at the surface of air and 
lossy metal or dielectric one has luvkz > 0, /ca, > 0; so 
that ImCp^x > 0 and < 0 and hence Eq. Q yields 
Pcirc,y < 0 (panel a). On the other hand, the field is 
actually incident upon the ultrathin wires from the side, 
rather than from the top face (panel b). As a result, 
the imaginary parts of the components of the field po¬ 
larization vector inside the wires are swapped (cp^x > 0 
and ImCp ^ < 0 ), and the sign of the circular polar¬ 
ization is reversed with respect to the continuous film. 
The crossover between the single wire and continuous 
film limit is illustrated by the polarization dependence on 
the lattice period shown in Fig. Qc). The polarization is 
positive for large periods in agreement with the geometric 
optics prediction (dotted line). For smaller periods the 
polarization is increased due to the collective effects and 





FIG. 4: Schematic illustration of the emergence of the circular 
polarization of difference signs for the p-polarized plane wave, 
incident upon (a) positively refracting medium and (b) wire 
array, (c) Dependence of averaged circular polarization inside 
the wires with D = 4.7 nm on the lattice period a. Dotted 
line shows the value for the a = 400 nm. 


it has a flat maximum at the spacings a ~ 10 4- 15 nm, 
close to the experimental data. For even smaller periods 
the polarization decreases and it changes the sign as a 
dense wire array approaches the continuous film limit. 

A good quantitative agreement between our theory and 
experiment is illustrated in Fig. [^e). In this figure we 
trace a ratio of the intensities of anti-Stokes and Stokes 
BLS lines as a function of the angle of incidence, both ex¬ 
perimental and theoretical values obtained using Eq. (|^ . 
In theoretical analysis we use the values of optical ellip- 
ticity obtained from our GST numerical simulations (see 
Eig. SI) while the same parameter for the SW is esti¬ 
mated using the analytical solution for Kittel modes on 
an individual ferromagnetic wire. To facilitate under¬ 
standing the angular evolution of the ellipticity of the 
optical polarization is given in the inset of Fig. [^e). 

In conclusion, by employing the Brillouin light scatter¬ 
ing tools we have observed the effect of anomalous po¬ 
larization conversion in arrays of ultrathin Go nanowires 
with the nanowire diameter below 5 nm manifested in the 
pronounced reversed Stokes/anti-Stokes scattering asym¬ 
metry. We have explained this effect by unexpected cir¬ 
cular polarization of light induced within the nanowires. 
In particular, the circular polarization has the opposite 
sign being much larger in the absolute value than that for 
the continuous films or thicker nanowires. This finding 
opens a great potential of seemingly simple nanowire ar¬ 
rays employed for manipulating light at the nanoscale. 
At the same time, our results suggest that the Bril¬ 
louin spectroscopy, traditionally employed as a probe of 
magnon states, is an extremely sensitive technique for 
studying polarization-resolved landscapes of local elec¬ 
tric fields in nanostructures. 
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Supplementary Information 

Here we present the details of the derivation of the re¬ 
sults for the spin mode frequency, that have been used 
in order to identify the dominantly scattering spin mode. 
In order to obtain the spin mode frequency as a func¬ 
tion of the wave vector q we average the Landau-Lifshitz 
equation of motion together with the Maxwell magnetic 
flux conservation over the cross-section of the cylinder, 
thus generalizing, for the cylindrical symmetry, the ap¬ 
proach, originally proposed for thin films by Stamps and 
Hillebrands [T] and revisited later in Ref. [5]. The ap¬ 
plicability of this technique, as mentioned before, is lim¬ 
ited to the case where the transversal distribution of the 
dynamic magnetization is close to uniform, which cor¬ 
responds to the Kittel mode we are interested in. This 
corresponds, in the Damon-Eshbach configuration, to the 
lowest fundamental mode, referred to in this article as 
“Kittel mode” 0. Below are sketched major steps in 
this calculation. A detailed account will be published 
elsewhere. 

A. The aim of this part is to derive the dynamic mag¬ 
netic field averaged over the section. In this section we 
use the coordinate system where the wire axis is directed 
along z. In the frame of the quasi static approximation, 
the dynamic magnetic field h is obtained as the gradient 
of a potential 

h = V$ . 

We assume a propagation along the nanowire axis z, 
^{x,y,z) = f{x,y)exp{iqz) . As the probed mode are 
quasiuniform across the section, we consider the dynamic 
magnetic field averaged over the section: 


where Kq is a modified Bessel function. 

B. In this part, the frequency from the averaged equa¬ 
tion of motion is derived assuming an effective anisotropy 
energy ttMJ -|- irMy — KM'^ /M'^, where K contains the 
magneto-crystalline contribution and the dipolar cou¬ 
pling. 

B.l First, we consider the case when the applied field is 
not saturating. The effective field reads Hy = H—2-KMy , 
Hz = 2KMz/M‘^. The equilibrium condition reads H = 
{2tt -|- 2K/M‘^)My. The averaged equations of motion 
yield 

i{n/j){ma;) = My[(hz)-2Aq‘^{mz)/M‘^+2K{mz)/M^] 
- Mz[{hy) - 2Aq^{my)/M'^] + {my)Hz - {mz)Hy . 


i(H/7)(TOy) = Mz[{K) - 2Aq^{m,;)IM'^] - {m,;)Hz , 


i(H/7)(m^) = -My[{h^) - 2Aq‘^{m^)/M'^] + {m^)Hy . 

Replacing the averaged field components by the expres¬ 
sions derived in the previous part, one obtains 
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B.2 Second, we consider the case the applied field is 
saturating. Using the method presented before, one ob¬ 
tains 


Next, we use the cylindrical coordinate system with ip 
being the azimuthal angle in the xy plane and p = 
\/a;2 -|- j/2 being the two-dimensional radius vector. Us¬ 
ing the Stokes theorem one can express the functions 
{d^/dx), (d^/dy), and iq($) via the values of f{p,p) 
evaluated at the nanowire surface where p = R. The 
function f{R, p) can be expanded as 

<f) = fo + fi cos p + f 2 smp . 

The coefficients /o, /i, /2 are derived from the boundary 
conditions: 

/o « (/), fi « -2ttR{ mx), /2 « -2ttR{ my) . 
Thus one obtains 

{hx) = -2Tr{mx), (hy) = -2Tr{my) , 


{h.) 
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C. In order to identify the origin of the dominant scat¬ 
tering mode we have measured the dependence of the 
Brillouin shift frequency D versus the applied magnetic 
field H; the results are presented in Fig. SI Interest¬ 


ingly, the fl{H) curve demonstrates a pronounced mini¬ 
mum: the frequency decreases for low values of H while 
displaying a clearly seen growth after passing the critical 
point. This characteristic softening is related to the sat¬ 
uration of the static magnetization. Red curve shows the 
calculation according to the theory above. The effective 
value of the wave vector q of the spin wave, contributing 
to the Brillouin scattering has been deduced by fitting the 
experimental results. According to these calculations the 
wavelength of the “scattering magnon” is 27r/(7 ~ 100 nm. 



















7 



* Electronic address: a.poddubny@phoi.ifmo.ru 

[1] R. L. Stamps and B. Hillebrands, Phys. Rev. B 44, 12417 
(1991) 

[2] F. Zighem, Y. Roussigne, S.-M. Cherif, and P. Moch, J 
Physics: Condensed Matter 19, 176220 (2007). 

[3] C.Kittel, Introduction to Solid State Phys. (Wiley, New 
York, 1996). 


FIG. SI: Dependence of the spin mode frequency on the ap¬ 
plied field strength. Squares are experimental results, the line 
has been calculated with the following set of parameters: M = 
1400 emu/cm^, A = 1.3x 10“®erg/cm, 7/(27r) = 3 GHz/kOe, 
R = 2.4 nm, K = —1.8 x 10® erg/cm®. 





